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combination of analytical and numerical results, the most entangled symmetric states for up to 12 
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Thomson's problem, and it is observed that, in general, they are different problems. 
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I. INTRODUCTION 

Being the fundamental resource in a wide range of situ- 
ations in quantum information processing, entanglement 
is considered as a 'standard currency' for quantum infor- 
mation tasks, and it is highly desirable to know which 
states of a given system exhibit a high or maximal amount 
of entanglement [Tj . When it comes to multipartite states 
this question becomes complicated. There are different 
types of entanglement [5] , alongside which there are many 
different ways to quantify entanglement, each of which 
may capture a different desirable quality of a state as a 
resource. 

In this work, the geometric measure of entanglement, a 
distance-like entanglement measure O H] , will be investi- 
gated to analyze maximally entangled multipartite states. 
There are several incentives to consider this particular 
measure. Firstly, it has a broad range of operational inter- 
pretations: for example, in local state discrimination |5j, 
additivity of channel capacities [5] and recently for the 
classification of states as resources for measurement-based 
quantum computation (MBQC)[7HSj. Another advantage 
of the geometric measure is that, while other known en- 
tanglement measures are notoriously difficult to compute 
from their variational definitions, the definition of the 
geometric measure allows for a comparatively easy calcu- 
lation. Furthermore, the geometric measure can be linked 
to other distance-like entanglement measures, such as the 
robustness of entanglement and the relative entropy of en- 
tanglement [H [ini E] . The function also has applications 
in signal processing, particularly in the fields of multi-way 
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data analysis, high order statistics and independent com- 
ponent analysis (IC A) , where it is known under the name 
rank one approximation to high order tensors {12H17j . 

We focus our attention on permutation-symmetric 
states ~ that is, states that are invariant when swapping 
any pair of particles. This class of states has been useful 
for different quantum information tasks (for example, in 
leader election \TB\). It includes the Greenberger-Horne- 
Zeilinger (GHZ) states [I9], W states and Dicke states 
|20j . and also occurs in a variety of situations in many- 
body physics. There has been lots of activity recently in 
implementing these states experimentally [2T1 [25] . Fur- 
thermore, the symmetric properties make them amenable 
to the analysis of entanglement properties [TOl . 

An important tool in this work will be the Majorana 
representation [25], a generalization of the Bloch sphere 
representation of single qubits, where a permutation-sym- 
metric state of n qubits is unambiguously mapped to n 
points on the surface of the unit sphere. Recently, the Ma- 
jorana representation has proved very useful in analyzing 
entanglement properties of symmetric states [23l |26l [27] . 
In particular, the geometric measure of entanglement has 
a natural interpretation, and the Majorana representation 
facilitates exploitation of further symmetries to character- 
ize entanglement |23j . For example, the two-qubit sym- 
metric Bell state = l/^/2 (|01) + |10)) is represented 
by an antipodal pair of points: the north pole |0) and the 
south pole Roughly speaking, symmetric states with 
a high degree of entanglement are represented by point 
distributions that are well spread out over the sphere. We 
will use this idea along with other symmetry arguments 
to look for the most entangled states. Along the way we 
will compare this problem to other optimization problems 
of point distributions on the sphere. 

The paper is organized as follows. In Section |ll| the 
definition and properties of the geometric measure of en- 
tanglement are briefly recapitulated, which is followed 
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by an introduction and discussion of symmetric states 



in Section III In Section IV the Majorana representa- 
tion of symmetric states is introduced. The problem of 
finding the maximally entangled state is phrased in this 
manner, and is compared to two other point distribution 
problems on S"^: Toth's problem and Thomson's prob- 
lem. In Section [Vj some theoretical results for symmetric 
states are derived with the help of the intuitive idea of the 
Majorana representation. The numerically determined 
maximally entangled symmetric states of up to 12 qubits 
are presented in Section [VI[ O ur results are discussed in 
Section [Vll[ and Section VIII| contains the conclusion. 



II. 



THE GEOMETRIC MEASURE OF 
ENTANGLEMENT 



The geometric measure of entanglement is a distance- 
like entanglement measure for pure multipartite states 
that assesses the entanglement of a state in terms of its 
remoteness from the set of separable states [H]. It is 
defined as the maximal overlap of a given pure state with 
all pure product states [H |30l El] and is also defined as 
the geodesic distance with respect to the Fubini-Study 
metric |32| . Here we present it in the inverse logarithmic 
form of the maximal overlap, which is more convenient in 
relation to other entanglement measures: 
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Eq is non-negative and zero iff \ip) is a product state. We 
denote a product state closest to l^/i) by |A^) € "Hsep, and 
it should be noted that a given \ip) can have more than one 
closest product state. Indeed, we will usually deal with 
entangled states that have several closest product states. 
Due to its compactness, the normalized, pure Hilbert 
space of a finite-dimensional system (e.g. n qudits) always 
contains at least one state j^*) with maximal entanglement, 
and to each such state relates at least one closest product 
state. The task of determining maximal entanglement 
can be formulated as a max-min problem, with the two 
extrema not necessarily being unambiguous: 
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It is often more convenient to define G{\tjj)) = 
max|;^) K'^IV')! : so that we obtain Eq = log2(l/G^). Be- 
cause of the monotonicity of this relationship, the task 
of finding the maximally entangled state is equivalent to 
solving the min-max problem 
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As mentioned in the Introduction, there are several advan- 
tages of this measure of entanglement. First of all, it has 
several operational interpretations. It has implications 
for channel capacity [B] and can also be used to give con- 
ditions as to when states are useful resources for MBQC 
[3j9]. If the entanglement of a set of resource states 
scales anything below logarithmically with the number of 
parties, it cannot be an efficient resource for determinis- 
tic universal MBQC |H]. On the other hand, somewhat 
surprisingly, if the entanglement is too large, it is also 
not a good resource for MBQC. If the geometric measure 
of entanglement of an n qubit system scales larger than 
n — S (where S is some constant), then such a computa- 
tion can be simulated efficiently computationally . Of 
course, we should also note that there are many other 
quantum information tasks that are not restricted by such 
requirements. For example, the n-qubit GHZ state can be 
considered the most non-local with respect to all possible 
two-output, two-setting Bell inequalities |33| , whereas the 
geometric measure is only ii^cdGHZ)) = 1, independent 
of n. The n-qubit W state, on the other hand, is the 
optimal state for leader election [TH| with entanglement 
Eg{\W)) = log2(n/(n - 1))""^ Indeed, for local state 
discrimination, the role of entanglement in blocking the 
ability to access information locally is strictly monotonic 
- the higher the geometric measure of entanglement, the 
harder it is to access information locally [5 . 

In addition, the geometric measure Eq has close links 
to other distance-like entanglement measures, namely the 
(global) robustness of entanglement R [34, and the relative 
entropy of entanglement i^R [5^ . Between these measures 
the inequalities Eq < Eb, < log2(l-l--R) hold for aU states 
in Uni [H], and they turn into equalities for stabilizer 
states (e.g. GHZ state), Dicke states (e.g. W state) and 
permutation-antisymmetric basis states [SJ [TOl [3S] . 

An upper bound for the entanglement of pure n qubit 
states is given in |36j as 



EQ{m<n~l 
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We can see that this allows for states to be more entangled 
than is useful, e.g. for MBQC. Indeed, although no states 
of more than two qubits reach this bound |36| , most states 
of n qubits have entanglement Eq > n — 21og2(n) — 3 
[7]. In the next section, we will see that symmetric states 
have generally lower entanglement. 

We can also make a general statement for positive 
states that will help us in calculating entanglement for 
this smaller class of states. For finite-dimensional systems, 
a general quantum state can be written in the form \ip) — 
ai\i) with an orthonormalized basis and complex 
coefficients G C. We will call \ip) a real state if - for a 
given basis {\i)} - all coefficients are real {ai G M), and 
likewise call {ip) a positive state if the coefficients are all 
positive {tti > 0). A computational basis is one made up 
of tensors of local bases. 

Lemma 1. Every state \^) of a finite-dimensional system 
that is positive with respect to some computational basis 
has at least one positive closest product state |A^). 
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Proof. Picking any computational basis in which the co- 
efficients of IV') are all positive, we denote the basis of 
subsystem j with {|ij)}, and can write the state as \ip) — 

SiajNi) ■ • • K")' with i = (zi, ...,?„) and aj> 0. A clos- 
est product state of can be written as |A^) = |fTj), 

where \aj) = }y'^_ (with hP^, e C) is the state of sub- 
system j. Now define a new product state with positive 
coefficients as |A;^) = (g)^. |ct^.), where = Y.^, Nj)- 

Because of |(^|A:^)| = Ya^YI, l&i-l > |E?«?n,^i = 
KV'IA^)!, the positive state |AJ^) is a closest product state 
of □ 

This lemma, which was also shown in ,37 , asserts that 
positive states have at least one positive closest product 
state, but there can nevertheless exist other non-positive 
closest product states. A statement analogous to Lemma 
[T] does not hold for real states, and it is easy to find 
examples of real states that have no real closest product 
state. 

From now on we will simply denote entanglement in- 
stead of referring to the geometric measure of entan- 
glement. It must be kept in mind, however, that the 
maximally entangled state of a multipartite system subtly 
depends on the chosen entanglement measure [38| . 



III. PERMUTATION SYMMETRIC STATES 

In general it is very difficult to find the closest prod- 
uct state of a given quantum state, due to the large 
amount of parameters in |A). The problem will be consid- 
erably simplified, however, when considering permutation- 
symmetric states. In experiments with many qubits, it 
is often not possible to access single qubits individually, 
necessitating a fully symmetrical treatment of the ini- 
tial state and the system dynamics [39]. The ground 
state of the Lipkin-Meshkov-Glick model was found to be 
permutation-invariant, and its entanglement was quanti- 
fied in term of the geometric measure and its distance- 
related cousins [4^. For these reasons it is worth analyz- 
ing various theoretical and experimental aspects of the 
entanglement of symmetric states, such as entanglement 
witnesses or experimental setups [4H I42j . 

The symmetric basis states of a system of n qubits are 
given by the Dicke states, the simultaneous eigenstates 
of the total angular momentum J and its z-component 
Jz [Sol EH |43] . They are mathematically expressed as 
the sum of all permutations of computational basis states 
with n — k qubits being |0) and k being |1). 



\Sn.k) = 



-1/2 



perm 



|0)|0).. -10)11)11) 



(5) 



with < fc < n, and where we omitted the tensor symbols 
that mediate between the n single qubit spaces. The Dicke 
states constitute an orthonormalized set of basis vectors 
for the symmetric Hilbert space Us- The notation \Sn,k) 



will sometimes be abbreviated as \Sk) when the number 
of qubits is clear. 

Recently, there has been a very active investigation 
into the conjecture that the closest product state of a 
symmetric state is symmetric itself [HUniUS]. A proof 
of this seemingly straightforward statement is far from 
trivial, and after some special cases were proofed [HI [3S] , 
Hiibener et al. were able to extend this result to the 
general case. They also showed that, for n > 3 qudits 
(general quantum d-level systems) the closest product 
state of a symmetric state is necessarily symmetric. This 
result greatly reduces the complexity of finding the closest 
product state and thus the entanglement of a symmetric 
state. 

A general pure symmetric state of n qubits is a linear 
combination of the n + l symmetric basis states, with the 
operational nature being that the state remains invariant 
under the permutation of any two of its subsystems. 

A closest product state of \Sn,k} is [10] 



|A) = 



n — k 



|o) 



ii; 



(6) 



i.e. a tensor product of n identical single qubit states. 
From this, the amount of entanglement is found to be 



EGi\Sn,k)) 




(7) 



This formula straightforwardly gives the maximally en- 
tangled Dicke state. For even n it is \Sn,n/2) and for odd 
n the two equivalent states |5'„_(„+i)/2) and |5'ra,(n-i)/2)- 
In general, however, the maximally entangled symmetric 
state of n qubits is a superposition of Dicke states. Never- 
theless, Equation ([t]) can be used as a lower bound to the 
maximal entanglement of symmetric states. This bound 
can be approximated by the Stirling formula for large n 
as Eq > log2 y/mr/2. 

An upper bound to the geometric measure for symmet- 
ric n qubit states can be easily found from the well-known 
decomposition of the identity on the symmetric subspace 
(denoted Isymm, see e.g. [55]), 



{\e){9\r-co{e) 
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^ Symm 
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where w denotes the uniform probability measure over 
the unit sphere S{H) on Hilbert space H. We can easily 

see that G(|V))' = niax<,e«,^p Tr(tj|V') (V-l) > 1/(^ + 1). 
Hence, for any symmetric state of n qubits, the geometric 
measure of entanglement is upper bounded by 



i^G(|^)s)<log2(n+l) 



(9) 



An alternative proof that has the benefit of being visually 
accessible is presented in Appendix [K\ 

The maximal symmetric entanglement for n qubits 
thus scales polylogarithmically between ©(log-y/n) and 
©(logn). To compare this with the general non-symme- 
tric case, consider the lower bound of the maximal n 
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qubit entanglement (n even) of i?G > ('^/2) [IZ]- Thus 
the maximal entanglement of general states scales much 
faster, namely linearly rather than logarithmically. As 
mentioned, for most states the entanglement is even higher 
and thus too entangled to be useful for MBQC. While 
the bounds for symmetric states mean that permutation- 
symmetric states are never too entangled to be useful 
for MBQC, unfortunately their scaling is also too low to 
be good universal deterministic resources [5]. They may 
nevertheless be candidates for approximate, stochastic 
MBQC [2]. Regardless of their use as resources for MBQC, 
the comparatively high entanglement of symmetric states 
still renders them formidable candidates for specific quan- 
tum computations or as resources for other tasks, such as 
the leader election problem [THj and LOCC discrimination 

We end this section by mentioning a simplification with 
respect to symmetric positive states. States that are 
symmetric as well as positive in some computational basis 
are labelled as positive symmetric. From the previous 
discussion it is clear that such states have a closest product 
state which is positive symmetric itself, a result first 
shown in [10] . It should be noted that, while each closest 
product state of a positive symmetric state is necessarily 
symmetric for n > 3 qudits, it need not be positive. We 
can formulate this as a statement akin to Lemma [H 

Lemma 2. Every symmetric state \'ijj)s of n qudits, which 
is positive in some computational basis, has at least one 
positive symmetric closest product state |A^)s. 



IV. MAJORANA REPRESENTATION OF 
SYMMETRIC STATES 



With the discussion of the geometric measure and sym- 
metric states behind us, we have gathered the prerequisites 
to introduce a crucial tool, the Majorana representation. 
It will help us to understand the amount of entanglement 
of symmetric states. 



A. Definition 

In classical physics, the angular momentum J of a 
system can be represented by a point on the surface of 
the 3D unit sphere S"^, which corresponds to the direction 
of J. No such simple representation is possible in quantum 
mechanics, but Majorana [28j pointed out that a pure 
state of spin-j can be uniquely represented by 2j not 
necessarily distinct points on S^. This is a generalization 
of the spin-1/2 (qubit) case, where the 2D Hilbert space 
is isomorphic to the unit vectors on the Bloch sphere. 

An equivalent representation also exists for permu- 
tation-symmetric states of n spin-1/2 particles [SHI fi5] . 
By means of this 'Majorana representation' any symmet- 
ric state of n qubits \ip)s can be uniquely composed from 
a sum over all permutations P : {1, . . . , 71} — > {1, . . . ,n} 



of n undistinguishable single qubit states {|0i), . . . , 

Ws^ -^^\(f>P(l))\(t>P(2))---\(t)p{n)) , (10) 
^ perm 

with 1 0,) = COS f 1 0) -t- e'"^' sin 1 1 1) , 



and K = nlY^ Y{{4'2\<Ppii)) ■ 

perm i—1 

The normalization factor K is in general different for 
different |V')s- By means of Equation (10), the multi- 
qubit state \Tp)s can be visualized by n unordered points 
(each of which has a Bloch vector pointing in its direction) 
on the surface of a sphere. We call these points the 
Majorana points (MP), and the sphere on which they lie 
the Majorana sphere. 

With Equation ( [l0| , the form of a symmetric state |'!/')s 
can be explicitly determined if the MPs are known. If the 
MPs of a given state |V')s = '^k=o'^k\Sk) are unknown, 
they can be determined by solving a system of 77, + 1 
equations. 

/n\ 

^ ( ) X] ^'P(l) ■ ' ' ^P{k)^P{k+l) ■ ■ ■ Cp(„) , (11) 
^ ^ perm 



with C; = COS ■ 



S, = e'"^' sin ■ 



The Majorana representation has been rediscovered 
several times, and has been put to many different uses 
across physics. In relation to the foundations of quantum 
mechanics, it has been used to find efficient proofs of 
the Kochen-Specker theorem [321 [50] and to study the 
'quantumness' of pure quantum states in several respects 
[5T| [52] , as well as the approach to classicality in terms 
of the discriminability of states [S3j. It has also been 
used to study Berry phases in high spin '54] and quantum 
chaos [551 156j . Within many-body physics it has been 
used for finding solutions to the Lipkin-Meshkov-Glick 
model [57", and for studying and identifying phases in 
spinor BEC [58 - 61^. It has also been used to look for 
optimal resources for reference frame alignment |62j and 
for phase estimation [63]. 

Recently, the Majorana representation has also become 
a useful tool in studying the entanglement of permutation- 
symmetric states. It has been used to search for and 
characterize different classes of entanglement [23l [26l [27] , 
which have interesting mirrors in the classification of 
phases in spinor condensates |23[ I59| . Of particular inter- 
est, in this work, is that it gives a natural visual inter- 
pretation of the geometric measure of entanglement |23J , 
and we will see how symmetries in the point distributions 
can be used to calculate the entanglement and assist in 
finding the most entangled states. 

The connection to entanglement can first be noticed 
by the fact that the point distribution is invariant under 
local unitary maps. Applying an arbitrary single-qubit 
unitary operation U to each of the n subsystems yields 
the LU map 



\(p)^ = [/(»••• (8)17 |V')s 



(12) 
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and from Equation ( 10 ) it follows that 



\f)s = I^P(l))I^P(2)) • • • I^^PH 



K 

perm 

with = u\(j)^) yi 



(13) 



In other words, the symmetric state |'0)s is mapped to 
another symmetric state \ip)s, and the MP distribution of 
\ip)s is obtained by a joint rotation of the MP distribution 
of |'i/')s on the Majorana sphere along a common axis. 
Therefore 1-0)8 and \(p)s have different MPs, but the same 
relative distribution of the MPs, and the entanglement 
remains unchanged. 

When it comes to the geometric measure of entangle- 
ment, we can be even more precise. For n > 3 qubits, 
every closest product state |A)s of a symmetric state 
is symmetric itself [25], so that one can write |A)s = \a) 
with a single qubit state \a), and visualize |A)s by the 
Bloch vector of |cr). In analogy to the Majorana points, 
we refer to \a) as a closest product point (CPP). 

For the calculation of the geometric measure of en- 
tanglement, the overlap with a symmetric product state 
|A) = |a)®" is 



(14) 



i=l 



The task of determining the CPP of a given symmet- 
ric state is thus equivalent to maximizing the absolute 
value of a product of scalar products. From a geometri- 
cal point of view, the are the angles between the 
two corresponding points on the Majorana sphere, and 
thus the determination of the CPP can be viewed as an 
optimization problem for a product of geometrical angles. 



B. Examples 

We will now demonstrate the Majorana representation 
for two and three qubit symmetric states. The case of 
two qubits is very simple, because any distribution of 
two points can be rotated on the Majorana sphere in 
a way that both MPs are positive, with |0i) = |0) and 
|(/,2) = cosf |0) -hsinf |1) for some 9 G [0,7r]. One CPP of 
this MP distribution is easily found to be jcr) = cos ||0) + 
sin 1 1 1). Figure 1 shows two examples for 6 = tt/2 and 
9 = TT, with the latter representing the Bell state = 
l/\/2(|01) + |10)). Due to the azimuthal symmetry of 
its MPs on the sphere, the CPPs form a continuous ring 
|ct) = 1/72(10) +e'^|l)), with if £ [0,27r) around the 
equator. The amount of entanglement is Eq{\iP'^)) = 1. 
For two qubits, the maximally entangled symmetric states 
are easily found to be those whose MPs lie diametrically 
opposite on the sphere. 

For three qubit states, the GHZ state and the W state, 
both of which are positive and symmetric, are considered 
as extremal among three qubit states [51]. The tripartite 




FIG. 1. The Majorana representations of two symmetric 
states of two qubits. MPs are shown as white circles and 
CPPs as dashed lines or crosses. Panel (a) depicts the state 
y^\00) + yi76(|01) + 1 10)). Its single CPP lies in the 
middle of the two MPs at an angle of 9 = 7r/4. Panel (b) 
shows the Bell state \ip+) = l/^/2(101) + |10)), whose CPPs 
form a continuous ring on the equatorial belt. 



GHZ state |GHZ) = 1/^2 (|000) + |111)) [TO] has the 
MPs 



6i) = ^(|0) + |l)) , 

i2) = ^(|0)+e'2-/3|l)) , 

^3>- 75(10) +e'^^/^|l)) . 



(15) 



Its two CPPs are easily calculated to be \ai) — |0) and 
\a2) = |1), yielding an entanglement of E'cdGHZ)) = 1. 
Figure [2ja) shows the distribution of the MPs and CPPs 
for the GHZ state. The three MPs form an equilateral 
triangle inside the equatorial belt, and the two CPPs lie 
at the north and the south pole, respectively. 

The W state |W) = 153,1) l/V^(|001)-h|010)-h|100)) 
is a Dicke state, and its MPs can be immediately accessed 
from its form as 



103 



> - \h) = |o) , 



ID 



(16) 



Generally, the definition ^ of the Dicke states \Sn,k) 
asserts that n — k MPs lie at the north pole and k at 




FIG. 2. The MPs and CPPs of the three qubit (a) GHZ 
state and (b) W state. The GHZ state has two discrete CPPs 
whereas for the W state the CPPs form a continuous ring due 
to the azimuthal symmetry. 
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the south pole. Equation (|6| yields \ai) = \/2/3 |0) + 
•\/l/3 |1) as a positive CPP of the W state, and from the 
azimuthal symmetry of the MP distribution it is clear 
that the set of all CPPs is formed by the ring of vectors 
\a) = ^/2/3\0) +e''fy/T/3\l), with if e [0,27r). Figure 
^h) shows the MPs and CPPs of |W). The amount of 
entanglement is Eq{\W)) — logj (9/4) 1.17, which is 
higher than that of the GHZ state. It was recently shown 
that, in terms of the geometric measure, the W state is 
the maximally entangled of all three qubit states 65J. 

It is insightful to examine how the CPPs and the entan- 
glement change when the MP distribution of the underly- 
ing state is modified. Starting out with three MPs lying 
on the north pole, two of the MPs are moved southwards 
on opposite sides (cf. Figure [s]), describing an isosceles 
triangle with the remaining MP on the north pole. Using 
the abbreviations ce = cos(0/2) and sg — sm{9/2), the 
MPs have the form 



= |0) , 
^2,3> =ce|0)±is(,|l) 



(17) 



with the parametrization 9 G [0, tt]. The form of the 



underlying quantum state follows from Equation ( 10 ) as 
3cj |000) -Hs^ (|011) + |101) + |110)) 



3si 



(18) 



This state is positive, so Lemma [2] asserts the existence 
of at least one positive CPP. With the ansatz jcr) — 



c^\0) +s^\l) for the CPP, the position of the CPP is 
found by calculating the absolute maximum of |(V'|cr)'*"^|. 
From this it is found that the parameter (p{9) of the CPP 
depends on the parameter 9 of the MPs as follows: 



l/{6sl - 3) 



(19) 



The permitted values of the left-hand side are [0,1], 
but the right-hand side lies outside this range for 9 < 
TT — arccos(l/5). For these values the CPP is fixed at 
|cr) = |0). Figure |3] shows how the CPP parameter ip{9) 
changes with 9. It is seen that from 9 = tt — arccos(l/5) 
onwards, the CPP abruptly leaves the north pole and 
moves towards the south pole along the prime meridian. 
From Equations ( 18 1 and ([T9| the amount of entangle- 



ment is easily calculated and is displayed in Figure [3] Eq 
is monotonously increasing [64] and reaches a saddle point 
at the GHZ state {9 = 27r/3). 



C. Entanglement and extremal point distributions 

The main point of interest in this paper is the study 
of maximally entangled symmetric states. For this the 
Majorana representation is extremely helpful, because 
it allows the optimization problem of maximizing the 
entanglement to be written in a simple form. With the 
help of Equation ( 14 ) , the min-max problem Q for finding 
the maximally entangled state can be reformulated as 




3 2 '5' 3 

MP position 9 

FIG. 3. Change of the entanglement and the location of the 
CPP when the MP distribution is modified. The position of the 
CPP does not change until the two moving MPs have reached 
a latitude slightly below the equator. From the distribution 
(c) onwards, the CPP rapidly moves southwards and reaches 
the equator at the GHZ state (d). After that, the location of 
the CPP and the entanglement changes only weakly until the 
W state (e) is reached. 




(20) 



This 'Majorana problem' bears all the properties of an 
optimization problem on the surface of a sphere in M.^. 
These kinds of problems deal with arrangements of a finite 
number of points on a sphere so that an extremal property 
is fulfilled |66j . Two well-known members, Toth's problem 
and Thomson's problem, have been extensively studied 
in the past. 

Toth's problem, also known as Fejes' problem and 
Tammes' problem, asks how n points have to be dis- 
tributed on the unit sphere so that the minimum distance 
of all pairs of points becomes maximal . This problem 
was first raised by the biologist Tammes in 1930 when try- 
ing to explain the observed distribution of pores on pollen 
grains [H7]. Recasting the n points as unit vectors e M'^, 
the following cost function needs to be maximized: 



/T6th(ri,r2, ■ • ■ ,rn) 



mm |r,; 



(21) 



The point configuration that solves this problem is called 
a spherical code or sphere packing |68j . The latter term 
refers to the equivalent problem of placing n identical 
spheres of maximal possible radius around a central unit 
sphere, touching the unit sphere at the points that solve 
Toth's problem. 

Thomson's problem, also known as the Coulomb 
problem, asks how n point charges can be distributed 
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on the surface of a sphere so that the potential energy 
is minimized. The charges interact with each other only 
through Coulomb's inverse square law. Devised by J. J. 
Thomson in 1904, this problem raises the question about 
the stable patterns of up to 100 electrons on a spherical 
surface [69j. Its cost function is given by the Coulomb 
energy and needs to be minimized. 



/Thomson(ri, r2, . . . , r„) 



i<j 



(22) 



The original motivation for Thomson's problem was to 
determine the stable electron distribution of atoms in 
the plum pudding model. While this model has been 
superseded by modern quantum theory, there is a wide 
array of novel applications for Thomson's problem or 
its generalization to other interaction potentials. Among 
these are multi-electron bubbles in liquid '*He [70j , surface 
ordering of liquid metal drops confined in Paul traps [71] , 
the shell structure of spherical viruses [751, 'colloidosomes' 
for encapsulating biochemically active substances ^73) . 
fuUerene patterns of carbon atoms [TJ] and the Abrikosov 
lattice of vortices in superconducting metal shells [75] . 

Exact solutions to Toth's problem are only known for 
riTo = 2—12, 24 points ^76], and in Thomson's problem 
for nTh = 2 — 8, 12 points [66ll76]. Despite the different 
definitions of the two problems, they share the same 
solutions for n = 2 — 6, 12 points [77 . Numerical solutions 
are, furthermore, known for a wide range of n in both 
problems [TSHST] . 



The solutions to n = 2, 3 are trivial and given by the 
dipole and equilateral triangle, respectively. For n — 
4,6,8, 12,20 the Platonic solids are natural candidates, 
but they are the actual solutions only for n = 4, 6, 12 
[5^ . For ri = 8, 20 the solutions are not Platonic solids 
and are different for the two problems. We will cover 
the solutions for n = 4 — 12 in more detail alongside the 
Majorana problem in Section \VJ\ 



On symmetry grounds, one could expect that the center 
of mass of the n points always coincides with the sphere's 
middle point. This is, however, not the the so- 

lution to Toth's problem for n = 7 [76] or the solution 
to Thomson's problem for n = 11 shows [7H1[75]. Fur- 
thermore, the solutions need not be unique. For Toth's 
problem, the first incident of this is n = 5 [53], and 
for Thomson's problem at n = 15 [75] and n = 16 [75] . 
These aspects show that it is, in general, hard to make 
statements about the form of the 'most spread out' point 
distributions on the sphere. The Majorana problem (20) 



is considered to be equally tricky, particularly with the 
normalization factor K depending on the MPs. Further- 
more, the MPs of the solution need not all be spread out 
far from each other, as demonstrated by the three qubit 
|W) state with its two coinciding MPs. 



V. STATES AND SYMMETRIES OF MP AND 
CPP DISTRIBUTIONS 

In this section, results for the interdependence between 
the form of n qubit symmetric states and their Majorana 
representation will be derived. More specifically, it will 
be examined what the distributions of MPs and CPPs 
look like for states whose coefficients are real, positive or 
vanishing. In some of these cases the MPs or CPPs have 
distinct patterns on the sphere, which can be described 
by symmetries. In this context, care has to be taken as 
to the meaning of the word 'symmetric'. Permutation- 
symmetric states were introduced in Section [TTl] and only 
these states can be represented by point distributions on 
the Majorana sphere. For some of these symmetric states, 
their MP distribution exhibits symmetry properties on the 
sphere. Examples of this can be found in Figure [2] where 
the GHZ state and W state have rotational symmetries 
around the Z-axis, as well as reflective symmetries along 
some planes. 

Let |V')s = J2k=oO,k\Sk) be a general symmetric state 
of n qubits. To understand the relationship between 
the state's coefficients and the Majorana representation, 
consider the effect of symmetric LUs. A symmetric LU 
acting on the Hilbert space of an n qubit system is defined 
as the n-fold tensor product of a single-qubit unitary 
operation: = U ^ ■ ■ ■ (E) U . This is precisely the LU 



map that was shown in Equation (12 1 and (13 1 to map 



every symmetric state to another symmetric state. 

Considering the Hilbert space of a single qubit, the 
rotation operator for Z-axis rotations of the qubit is 



1 

e'*^ 



(23) 



and the rotation operator for F-axis rotations of the qubit 
is 



Ry{9) = 



(24) 



i?2 changes the relative phase, but not the absolute value 
of the qubit's coefficients. Conversely, Ry changes the ab- 
solute value, but not the relative phase of the coefficients. 
From Equation ( 10 1 it is easily seen that Rz and R 



this behavior on to the symmetric LUs R^ 



R,] 



For example, the effect of i?| on 



RmWs 



fc=0 



Sk) 



y pass 
and 



(25) 



From this it is easy to determine the conditions for the 
MPs of \'ip)s having a rotational symmetry around the 
Z-axis, i.e. i?|(0) |V')s = \'4')s (up to a global phase) for 
some 6 < 2tt. From Equation (25 1 it is clear that the 



possible rotational angles (up to multiples) are restricted 
to 9 = 27r/m, with m G N, 1 < m < n. The necessary 
and sufficient conditions are: 
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Lemma 3. The MP distribution of a symmetric n qubit 
state |^)s is rotationally symmetric around the Z-axis 
with rotational angle 6 = 27r/m (1 < m < n) iff 

kj\ Uk, ^OA akj 7^ 0} : {h - kj) mod m^O (26) 



Proof. Equation (26) is equivalent to: 3 Z G Z : 
V{fc| ak ^ 0} : k mod m = I. From this it fol- 
lows that i?|(27r/m) |'0)s = exp(i27rfc/m)|S'fe) = 



exp(i27r//m)|iS'fc) — e |V')si with 5 — 2-1:1/7 



□ 



In other words, a sufficient number of coefficients need 
to vanish, and the spacings between the remaining coeffi- 
cients must be multiples of m. For example, a symmetric 
state of the form |-0)s = aslSa) -I- arjiSy) -I- oisjS'is) is ro- 
tationally symmetric with 9 — 7r/2, because the spacings 
between non-vanishing coefficients are multiples of 4. 

Lemma 4. Every maximally entangled symmetric state 
\^) s of n qubits has at least two different CPPs. 

Proof. The cases n = 2, 3 are trivial, because their maxi- 
mally entangled states (Bell states and W state, respec- 
tively) have an infinite number of CPPs. For n > 3, we 
consider a symmetric state \tjj) with only one CPP \a) 
and show that \ip) cannot be maximally entangled. 

Because of the LU invariance on the Majorana sphere 
(cf. Equation (13)), we can take the CPP to be the 
north pole, i.e. \a) = |0). Denoting a single qubit with 
\Ld) = C6i|0) -f e"^S6(|l), the smooth and continuous overlap 
function g{\uj)) = |(V'|a;)'*"| then has its absolute max- 
imum at |cj) = |0). For any other local maximum |a;'), 
the value of g{\ijj')) is smaller than g{\0)) and therefore 
an infinitesimal change in the MPs of \tp) cannot lead to 
a CPP outside a small neighborhood of the north pole. 

We will now present an explicit variation of \ip) that 
increases the entanglement. 1-0) = '^k=oO'k\Sk) has 
complex coefficients that fulfil (V'l^) — 1, as well as 
ao > and Oi = [S3]. Define the variation as = 
(ao-e)|5o)+I]fc=2 a/c|S'fe)-l-(a„+eao/a* )|S'„), with e < 1. 

This state fulfils the requirement Itpa) \tp), and is nor- 
malized: (V'fl'i/'e) = l-l-C'(e^). We now investigate the val- 
ues of ged^)) = |(V'e|w)®"| around the north pole. In this 
area = (1 - {0y8) + O{9^))\O)+e'^{{0/2)-O{e^))\l), 
hence g^u;)) = |(1 - (0V8))"(«o - e)\ + 0{e\9^) = 
|ao — e| +0(6^,0^) < |ao| — g(|0)) for smaU, but nonzero 
e and 9. Therefore the absolute maximum of is smaller 
than that of g, and is more entangled than \-4>). □ 



A. Real symmetric states 

For symmetric states with real coefficients, the following 
lemma asserts a reflection symmetry of the MPs and CPPs 
with respect to the X-Z-plane that cuts the Majorana 
sphere in half. In mathematical terms, the MPs and 
CPPs exhibit a reflection symmetry with respect to the 
X-Z-plane iff for each MP = ce|0) -I- e^'^sell) the 
complex conjugate point \(f)i)* = Cg\Q) + e~"''se|l) is also 
a MP, and the same holds for CPPs too. 



Lemma 5. Let s be a symmetric state of n qubits. 
lip) s is real iff all its MPs are reflective symmetric with 
respect to the X -Z -plane of the Majorana sphere. 

Proof. (^) Let |i/')s be a real state. Then |V')s = Wt, 
and since Majorana representations are unique, |'0)s has 
the same MPs as |V')s- Therefore the complex conjugate 
of each non-real MP \(pi) is also a MP. 
(-4=) Let the MPs of |V')s be symmetric with re- 
spect to the X-Z-plane. Then for every nonreal MP 
its complex conjugate \4>i)* is also a MP. Because 
{\(j)i)\(f>i)* + \4>i)*\4'i)) is real, it becomes clear, from the 
permutation over all MPs in Equation (10), that the 
overall state \ip)s is real, too. □ 

The reflective symmetry of the MPs naturally leads to 
the same symmetry for the CPPs. 

Corollary 1. Let \tjj)s be a symmetric state of n qubits. 
Lf \ip)s is real, then all its CPPs are reflective symmetric 
with respect to the X -Z -plane of the Majorana sphere. 

Proof. Lemma [5] asserts that for every MP of |V')si the 
complex conjugate \4>i)* is also a MP. By considering the 
complex conjugate of the optimization problem (20), it 



becomes clear that for any CPP \a) the complex conjugate 
I cr)* is also a CPP. □ 



B. Positive symmetric states 

For symmetric states with positive coefficients, strong 
results can be obtained with regard to the number and 
locations of the CPPs. In particular, for non-Dicke states 
it is shown that there are at most 2n — 4 CPPs and that 
non-positive CPPs can only exist if the MP distribution 
has a rotational symmetry around the Z-axis. Further- 
more, the CPPs can only lie at specified azimuthal angles 
on the sphere, namely those that are 'projected' from the 
meridian of positive Bloch vectors by means of the Z-axis 
rotational symmetry (see, e.g., the positive seven qubit 
state shown in Figure [t]). 

Dicke states constitute a special case due to their con- 
tinuous azimuthal symmetry. The two Dicke states |5o) 
and \Sn) are product states, with all their MPs and CPPs 
lying on the north and the south pole, respectively. For 
any other Dicke state \Sk) the MPs are shared between 
the two poles, and the CPPs form a continuous horizontal 
ring with inclination 9 — 2 arccos \Jn — k/n. 



Lemma 6. Let \^l))s be a positive symmetric state of n 
qubits, excluding the Dicke states. 

(a) If s is rotationally symmetric around the Z-axis 
with minimal rotational angle 27r/m, then all its 
CPPs \cT{9,ip)) = ce|0) -I- e^'fsg\l) are restricted to 
the m azimuthal angles given by (f — (fr — 2'Kr/m 
with r (zli. Furthermore, if \a{9,ipr)) is a CPP for 
some r, then it is also a CPP for all other values 
ofr. 
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(b) If \ip)s is not rotationally symmetric around the 
Z-axis, then all its CPPs are positive. 

Proof. The proof runs similar to the one of Lemma [l| 
where the existence of at least one positive CPP is es- 
tablished. We use the notations |^)s = '^f.ak\Sk) with 
ak > 0, and |A) = |cr)®". 

Case (a): Consider a non-positive CPP \(7) — cg|0) + 
e^^Sejl) with if = 2'Kr/m and r G K, and define |cr) — 

ikip 



ce|0)+se|l). ThenKAIV'^sl = I Efc e'^^'^OfeC^ 
J2k o-kCg~'^Sgy/n/k — |(A|V')s|- If this inequality is strict, 
then |cr) is not a CPP. This would be a contradiction, so 
it must be an equality. Thus, for any two indices ki and 
kj of non-vanishing coefficients Ofc. and ak , the following 




must hold: e 

kjr mod m = fc. 



= e^^"^. This can be reformulated as 
or equivalently 



kj ) r mod m 







(27) 



Because ip = 2it /m is the minimal rotational angle, to is 
the largest integer that satisfies Equation (26), and thus 



there exist ki and kj with Ofc . , ^ s.t. ki — kj ^ 



From this and from Equation (27), it follows that r G Z 



Therefore |cr(0, ipr)) is a CPP if and only if r is an integer. 

Case (b): Considering a CPP |cr) = ce|0) + e'''se|l) 
with p = 27rr and r € M, we need to show that r e Z. 
Defining \a) = Cg\0) + Sg\l), and using the same hue of 



kip 



must 



argumentation as above, the equation e' 
hold for any pair of non- vanishing a^- and a^^ . This is 
equivalent to 



{ki — kj) r mod 1 = 



(28) 



or {ki—kj)r g Z, in particular r £ Q. If there exist indices 



ki and kj of non-vanishing coefficients s.t. 



ki kj — 1 , 



then r € Z, as desired. Otherwise consider r = x/y with 
coprime a;, j/ € N, a; < y. Because |V')s is not rotationally 
symmetric, the negation of Lemma [3] yields that, for any 
two ki and kj (a^^, a^. ^ 0) with ki — kj ~ a > 1, there 
must exist a different pair kp and kq (ofe^, Ok^ 0) with 
kp~kq = j3 > \ s.t. a is not a multiple of /3 and vice versa. 



From r = x/y and Equation (28), it follows that y = a 
as well as y = f3. This is a contradiction, so r e Z. □ 

With this result about the confinement of the CPPs 
to certain azimuthal angles, it is possible to derive upper 
bounds on the number of CPPs. 

Theorem 1. The Majorana representation of every posi- 
tive symmetric state \ip)s of n qubits, excluding the Dicke 
states, belongs to one of the following three mutually ex- 
clusive classes. 

(a) \'ip)s is rotationally symmetric around the Z-axis, 
with only the two poles as possible CPPs. 

(b) \'ip)s is rotationally symmetric around the Z-axis, 
with at least one CPP being non-positive. 

(c) s is not rotationally symmetric around the Z-axis, 
and all CPPs are positive. 



Regarding the CPPs of states from class (b) and (c), the 
following assertions can be made for n > 3; 

(b ) If both poles are occupied by at least one MP each, 
then there are at most 2n — 4 CPPs, else there are 
at most n CPPs. 

(c) There are at most [^^] CPPs J^J. 

Proof. Starting with the first part of the theorem, case (c) 
has already been shown in Lemma [6j so consider states 
\'ip)s with a rotational symmetry ip = 2T:/m around the 
Z-axis. If all CPPs are either |0) or then we have 
case (a), otherwise there is at least one CPP \a) which 
does not lie on a pole. If this |cr) is non-positive, then we 
have case (b), and if |cr) is positive, then Lemma |6] states 
the existence of another, non-positive CPP, thus again 
resulting in case (b). 

The proof of the second part of the theorem is a bit 
involved and can be found in Appendix |B] □ 



VI. NUMERICAL SOLUTIONS FOR UP TO 
TWELVE QUBITS 

In this section, we present the numerically determined 
solutions of the Majorana problem for up to 12 qubits. In 
order to find these, the results from the previous sections 
were extremely helpful. This is because an exhaustive 
search over the set of all symmetric states quickly becomes 
unfeasible, even for low numbers of qubits and because 
the min-max problem ( 20 ) is too complex to allow for 



straightforward analytic solutions. 

Among the results particularly helpful for our search 
are Lemma |6] and Theorem [l] For positive states they 
strongly restrict the possible locations of CPPs, and thus 
greatly simplify the calculation of the entanglement. It 
then suffices to determine only the positive CPPs be- 
cause all other CPPs automatically follow from the Z-axis 
symmetry (if any exists). We will see that this result is 
especially powerful for the Platonic solids in the cases 
n = 4, 6 where the location of the CPPs can be immedi- 
ately determined from this argument alone, without the 
need to do any calculations. 

From the definition ^ of Eq, it is clear that the exact 
amount of entanglement of a given state (or its corre- 
sponding MP distribution) is automatically known once 
the location of at least one CPP is known. A numerical 
search over the set of positive symmetric states often de- 
tects the maximally entangled state to be of a particularly 
simple form, enabling us to express the exact positions of 
its MPs and CPPs analytically. In some cases, however, 
no analytical expressions were found for the positions of 
the CPPs and/or MPs. In these cases the exact amount 
of entanglement remains unknown, although it can be 
numerically approximated with high precision. 

In this way we can be quite confident of finding the 
maximally entangled positive symmetric state. In the 
general symmetric case we do not have as many tools, so 
the search is over a far bigger set of possible states, and we 
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can be less confident in our results. We therefore focus our 
search to sets of states that promised high entanglement. 
Such states include those with highly spread out MP 
distributions, particularly the solutions of the classical 
optimization problems. From these results we can propose 
candidates for maximal symmetric entanglement. 

For two and three qubits, the maximally entangled 
symmetric states are known and were discussed in Section 
IV so we start with four qubits. A table summarizing the 
amount of entanglement of the numerically determined 
maximally entangled positive symmetric as well as the 
entanglement of the candidates for the general symmetric 
case can be found in Appendix [C| 



A. Four qubits 

For four points, both Toth's and Thomson's problem 
are solved by the vertices of the regular tetrahedron 
The numerical search for the maximally entangled sym- 
metric state returns the Platonic solid too. Recast as 
MPs, the vertices are 



I'/'i) - |0) , 

l'/'4 



- J5lo> + Viii) 



V3I 

73 



,i27r/3 



+ C 



§11) 
i4V3 7||l) 



(29) 



The symmetric state constructed from these MPs by 
means of Equation ( 10 ) shall be referred to as the 'tetra- 
hedron state'. 



Its form is 1*4) = Vl/3 \So) + jS'a), 
and its MP distribution is shown in Figure [4] Because the 
state is positive and has a rotational symmetry around 
the Z-axis, Lemma [6] restricts the possible CPP locations 
to the three half-circles \<j{9, ip)) with tp = 0, 27r/3, Att/S. 
From the symmetry of the Platonic solid it is clear that 
the MP distribution of Figure |4] can be rotated s.t. \(f>2), 
14)3) or |(j!)4) is moved to the north pole, with the actual 
distribution (and thus |^'4)) remaining unchanged. Each 
of these rotations, however, gives rise to new restrictions 
on the location of the CPPs mediated by Lemma |6] The 
intersections of all these restrictions are the four points 
where the MPs lie. This yields the result that 1^*4) has 




four CPPs, with their Bloch vectors being the same as 
those in Equation ( 29 ) . Fron 
glement follows as i?G(|^'4)) = log2 - 



those in Equation (|29|) . From this the amount of entan- 

,3« 1.59. 



B. Five qubits 

For five points, the solution to Thomson's problem is 
given by three of the charges lying on the vertices of an 
equatorial triangle and the other two lying at the poles 
Emu]- This is also a solution to Toth's problem, but it 
is not unique [83l[87j. The corresponding quantum state, 
the 'trigonal bipyramid state', is shown in Figure [5](a) . 
This state has the form \ip5) = l/V2{\Si) + {Si)), and 
a simple calculation yields that it has three CPPs that 
coincide with the equatorial MPs, giving an entanglement 
ofi;G(|V^5)) = log2(16/5)«1.68. 




FIG. 5. The distribution (a) shows the 'trigonal bipyramid 
state', but the conjectured solution of the Majorana problem 
is the 'square pyramid state', shown in (b). 

However, a numerical search among all positive sym- 
metric states yields states with higher entanglement. Our 
numerics indicate that the maximally entangled state is 
the 'square pyramid state' shown in Figure [Sjb). This 
state has five CPPs, one on the north pole and the other 
four lying in a horizontal plane slightly below the plane 
with the MPs. The form of this state is 



1*5) = 



\So) 



A 



\S4 



Its MPs are 



I0i) = |0) , 



102,3,4,5) =a|0)+c'Vl-a2|l) _ 

with K = |, ^, ^, and the CPPs are 

Wi) - |o) , 

k2,3,4,5) -xlO) +fcv/l-a;2|l) , 



(30) 



(31) 



(32) 



FIG. 4. MPs and CPPs of the 'tetrahedron state'. 



with fc = 1, i, — 1, — i. The exact values can be determined 
analytically by solving quartic equations. The x € (0, 1) 
of Equation (|32| is given by the real root of 4a;^ -I- ix^ + 
4a;2 — a; — 1 = 0, and this can be used to calculate A ~ 
(1 - x^)/{V5x{l - x^ f). With the substitution a = 
c? E (0, 1) the value of a is given by the real root of 
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(5A2-l)a4 + 4a3-6a2- 
of these quantities are: 



- 4a — 1 = 0. Approximate values 



0.46657 



A « 1.53154 



0.59229 



(33) 



The entanglement is iJcd^'s)) = log2(l + ~ 1-74, 
which is considerably higher than that of Eq^i/j^)). We 
remark that the 'center of mass' of the five MPs of |\E'5) 
does not coincide with the sphere's origin, thus ruling out 
the corresponding spin-5/2 state as an anticoherent spin 
state, as defined in [5T] . 



C. Six qubits 

The regular octahedron, a Platonic solid, is the unique 
solution to Toth's and Thomson's problem for six points. 
The corresponding 'octahedron state' was numerically 
confirmed to solve the Majorana problem for six qubits. 











- \4>e)-'^ 










I04> 












FIG. 6. Two possible orientations of the 'octahedron state'. 

The straightforward orientation shown in Figure [6]^a) 
has the form 1^'^,) = l/^/2(\S^) - jS's)), and its MPs are 



\%) = l/^/2(|5l) 

l0l) = |O), \cl)2 
103,4,5,6) = ^(|( 



(34) 



with fc = 1, i, —1, — i. I^'g) can be turned into the positive 
state 1*6) = l/V2{\Si) + l^s)) by means of an i?|(7r/4) 
rotation. The CPPs can be obtained from this state 
in the same way as for the tetrahedron state. Being a 
Platonic solid, the MP distribution of Figure [6](b) is left 
invariant under a finite subgroup of rotation operations 
on the sphere. From Lemma |6| the intersection of the 
permissible locations of the CPPs is found to be the eight 
points lying at the center of each face of the octahedron, 
forming a cube inside the Majorana sphere. 



Fl,2,3,4/ 
k5,6,7,8) 



10^ 

2^3 



-■/^|o) 



2^3 1-^/ 



V3+1 
2\/3 



(35) 



1) 



with k = l,i, — 1,— i. In contrast to the tetrahedron 
state, where the MPs and CPPs overlap, the CPPs of the 
octahedron state lie as far away from the MPs as possible. 
This is plausible, because in the case of the octahedron 



Equation ( 20 1 is zero at the location of any MP, due to the 
MPs forming diametrically opposite pairs. The amount 
of entanglement is Scd^'e)) = log2(9/2) « 2.17. 



D. Seven qubits 

For seven points, the solutions to the two classical 
problems become fundamentally different for the first time. 
Toth's problem is solved by two triangles asymmetrically 
positioned about the equator and the remaining point 
at the north pole [76], or (1-3-3) in the Foppl notation 
[55] . Thomson's problem is solved by the vertices of a 
pentagonal dipyramid [751 HH US] , where five points lie 
on an equatorial pentagon and the other two on the poles. 
The latter is also numerically found to be the solution to 
the Majorana problem. 




FIG. 7. MPs and CPPs of the 'pentagonal dipyramid state'. 
The ten CPPs are equidistantly located on two circles above 
and below the equator. 

The 'pentagonal dipyramid state', shown in Figure [7| 
has the form l^-y) = 1/^2(1 5i) + 15*6)), and its MPs are 



|0i) = |O), 102) = |1) , 
103,4.5,6,7) = ;)3(|0)+c''^|l)) 



(36) 



with K = 0, ^, ^, ^, ^. The CPPs of this positive state 
can be determined analytically by choosing a suitable 
parametrization. With x :— cos^ 9 the position of |cti) = 
ce|0) -I- Se|l) and |(T2) — Sg\0) -I- c^jl) is determined by 
the real root of the cubic equation 49a;'^ -I- 165x^ — 205x -|- 
55 = in the interval [0, The approximate amount of 
entanglement is £;g(|*7)) ~ 2.299. 



E. Eight qubits 

For eight points, Toth's problem is solved by the cu- 
bic antiprism, a cube with one face rotated around 45° 
and where the distances between neighboring vertices 
are all the same. The solution to Thomson's problem 
is obtained by stretching this cubic antiprism along the 
Z-axis, thereby introducing two different nearest-neighbor 
distances and further lowering symmetry [76l [78l [86] . One 
would expect that a similar configuration solves the Ma- 
jorana problem too, but, surprisingly, this is not the case. 
The 'asymmetric pentagonal dipyramid' shown in Figure 
|8]jb) is numerically found to have the highest amount of en- 
tanglement. An analytic form of this positive state is not 
known, but it can be numerically approximated to a very 
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high precision. The state is l^s) » 0.672|S'i) + 0.741|S'6), 
and its entanglement is iScd^s)) ~ 2.45. For comparison, 
the regular cube yields -E'cdV'cubc)) — log2(24/5) « 2.26. 
Interestingly, the positive state j^'g) has a higher amount 
of entanglement than any state of the antiprism form, all 
of which are non-positive. Furthermore, two of the MPs 
of l^fg) coincide, akin to the W state of three qubits. 



inclination 6 of the CPPs in the northern hemisphere (or 
TT ~ 9 for the CPPs in the southern hemisphere) follows 
from the real root of 81a;^ + 385x^ - 245a; + 35 = in the 
interval [0, 0.3]. The approximate amount of entanglement 
is ^g(|*9)) « 2.554. 














) — 'Tl/T. 







FIG. 8. The maximally entangled antiprism state is shown in 
(a), while the 'asymmetric pentagonal dipyramid state' in (b) 
is conjectured to be the maximally entangled state. 

I^'s) is positive, with two MPs lying on the north pole, 
one on the south pole and the other five on a circle below 
the equator. The exact inclination of this circle as well 
as the inclination of the two circles with the CPPs is not 
known, but can be approximated numerically. 



F. Nine qubits 

For nine points, the solutions to Toth's and Thom- 
son's problems are slightly different manifestations of 
the same geometric form (3-3-3) of neighboring triangles 
being asymmetrically positioned. This is also known 
as a triaugmented triangular prism. In contrast to 
this, the Majorana problem is numerically solved by 
1*9) = 1/72(152) + l^y)), shown in Figure § This is 
again a positive state with a rotational Z-axis symmetry 
and with coinciding MPs. 

The CPPs can be found analytically in the same way 
as for seven qubits. With the substitution x := cos^ 6, the 




G. Ten qubits 

The solution to Toth's problem is an arrangement of 
the form (2-2-4-2), while Thomson's problem is solved by 
the gyroelongated square bipyramid, a deltahedron that 
arises from a cubic antiprism by placing square pyramids 
on each of the two square surfaces. 

The ten-qubit case is distinct in two respects. It is 
the first case where the numerically determined positive 
solution is not rotationally symmetric around any axis. 
Furthermore, we found non-positive states with higher 
entanglement than the conjectured solution for positive 
states. 

A numerical search returns a state of the form |'3'io) = 
alSo) + PlSi) -1-7159) as the positive state with the highest 
entanglement, namely £'g(|^'io)) ~ 2.6798. From Lemma 
[3] it is clear that this state is not rotationally symmetric 
around the Z-axis. The MP distribution is shown in Fig- 
ure [TojJ a). The state has only three CPPs, which are all 
positive (cf. Theorem [T]), but there are six other local 
maxima of g{a) with values close to the CPPs. Their posi- 
tions are shown by dashed crosses in Figure 10 'a) . While 
the total MP distribution is not rotationally symmetric 
around the Z-axis, one would expect from the numerical 
results that the MPs form two horizontal planes, one with 
five MPs and another with four MPs, with equidistantly 
spread out MPs. However, this is not the the 
locations of the MPs deviate by small, but significant, 
amounts from this simple form. 




FIG. 10. The numerically determined maximally entangled 
positive state is shown in (a). A similarly highly entangled 
positive state with a rotational symmetry is shown in (b) . The 
candidate for the general case is shown in (c). 




FIG. 9. The 'pentagonal dipyramid state' with both of the 
poles occupied by two MPs. 



Interestingly, there is a fully rotationally symmetric 
positive state that comes very close to |^io) in terms 
of entanglement. Its straightforward form is IV'io) — 
1/72(152) 1 58)), as displayed in Figure floFb). The 12 



CPPs of this state are easily found as the solutions of a 
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quadratic equation. The two positive CPPs are 



I. Twelve qubits 



1 

V3+\/3 



|0) 
|0> 



(37) 



\/3-\/3 



and the entanglement is Eq{\^iq)) = log2(32/5) w 2.6781. 
This is less than 0.1% difference from |4'io)- 

The solution to Thomson's problem, recast as a quan- 



tum state of the form 



10/ 



a\S{)+p\S^)-a\S^), 



not positive and has an entanglement of Eq sa 2.7316. 
Prom numerics one can see that the entanglement of this 
state can be further increased by slightly modifying the 
coefficients, arriving at a state with eight CPPs and an 
entanglement of iScd^'io)) ~ 2.7374. The state is shown 
in Figure [To|^c) , and we propose it as a candidate for the 
maximally entangled symmetric state of ten qubits. 



For 12 points, both of the classical problems are solved 
by the icosahedron, a Platonic solid. Because the icosa- 
hedron cannot be cast as a positive state, the numerical 
search for positive states yields a different state of the 
form \%^) = a|S'i)+/3|S'6)+a|S'ii). From Figure [l2];a) it 
can be seen that this state can be thought of as an icosa- 
hedron with one circle of MPs rotated by 36° so that it is 
aligned with the MPs of the other circle. There are three 
circles of CPPs with five in each circle. One of these planes 
coincides with the equator, so that \a) = 1/a/2(|0) + |1)) 
is a trivial CPP. Nevertheless, the exact location of some 
of the MPs and CPPs are unknown. The approximate 
entanglement is £;G(|*i2)) ~ 2.99. 





H. Eleven qubits 

The solution to Toth's problem is a pentagonal an- 
tiprism with a pentagonal pyramid on one of the two 
pentagonal surfaces, or (1-5-5). The solution to Thom- 
son's problem is of the form (1-2-4-2-2). Analogous to 
the ten-qubit case, the numerically found positive state of 
11 qubits with maximal entanglement is not rotationally 
symmetric. The state, shown in Figure 11 |a), is of the 
form 1^11) — a\Si) + (3\S^) +7|5io), and its entanglement 
is £'g(|*ii)) ~ 2.77. The state has only two CPPs, but 
there exist seven more local maxima with values close to 
the CPPs. 




FIG. 11. The conjectured maximally entangled positive state 
of 11 qubits is shown in (a), while the candidate for the general 
case is shown in (b). 

The solution to Toth's problem, which is of the form 
IV'n) = o\Sq) + I3\S^) - 7|S'io), yields very low entangle- 
ment, but by modifying the coefficients of this non-positive 
state one can find a state which is even more en- 

tangled than 1^11). As shown in Figure [Tl|(b), the state 
is rotationally symmetric around the Z-axis and has 11 
CPPs. The entanglement is £:G(|*ii)) ~ 2.83, making 
the state the potentially maximally entangled state of 11 
qubits. 




FIG. 12. An orientation of the maximally entangled positive 
symmetric state of 12 qubits is shown in (a). The icosahe- 
dron state, shown in (b), is considered to be the maximally 
entangled of all symmetric 12 qubit states. 

Due to the high symmetry present in Platonic solids, 
the 'icosahedron state' is a strong candidate for max- 
imal symmetric entanglement of twelve qubit states. 
The state can be cast with real coefficients [^'12) = 



^l-^i) - ^\Sq) - ^\Sii), and its MPs can be eas- 
ily derived from the known angles and distances in the 
icosahedron. 



l'/'l) = |0), \<t>12)^\l) 
l<^2,3,4,5,6 

|</>7,8,9,10,li; = 



5+\/5 



(38) 



with K = 0, From numerics and from the 

Z-axis rotational symmetry, it is evident that there are 
20 CPPs, one at the center of each face of the icosahe- 
dron. Equivalent to the six-qubit case, the MPs appear 
as diametrically opposite pairs, forcing the CPPs to be 
as remote from the MPs as possible. The CPPs are 



ki6,.. 



.,5) =a+|0)+e*('=+'^/5) a_|l) 
,10) -b+|0)+c'(-+-/5)b_|l) 
,15) =b_|0)+e"'b+|l) , 
,20) =a_|0)+e'-a+|l) , 



(39) 
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with K = 



27r 
5 ' 



4tv 
5 ' 



5 ' 



5 



and 



a± 



1 1 /5 + 2V^ 



b+ = 



(40) 



\ 2 2 V 15 



With the knowledge of the exact positions of the MPs 
and CPPs, the entanglement of the icosahedron state 
can be calculated as £^0(1*12)) = log2(243/28) « 3.1175. 
Figure |13| shows a spherical plot of the overlap function 
9i'^) — l(^i2|o')'^^^| from the same viewpoint as in Figure 
[l2|[b). Due to their diametrically opposite pairs, the MPs 
coincide with the zeros in this plot. The CPPs can be 
identified as the maxima of (?(cr). 




FIG. 13. (color online) A spherical plot of the overlap function 
g(a) = |(*i2|o-)®^^| for the icosahedron state 1*12). 



VII. DISCUSSION 

The MP distribution of highly entangled states can 
be explained with the overlap function g{a) = |(V'|o')'*"| 
seen in Figure [73} Appendix states that the integra- 
tion volume of gicr)^ over the sphere is the same for all 
symmetric states. Therefore a bunching of the MPs in 
a small area of the sphere would lead to high values of 
g{<jY in that area, thus lowering the entanglement. This 
explains the tendency of MPs to spread out as far as 
possible, as it is seen for the classical problems. How- 
ever, there also exist highly entangled states where two 
or more MPs coincide (as seen for n — 3,8,9). This 
is intriguing because such configurations are the least 
optimal ones for classical problems. Again, this can be 
explained with the constant integration volume of g{(7)^. 
Because of 5(0-)^ cx J|- |((^i|(T)p, the zeros of g{a)'^ are 
the diametrically opposite points (antipodes) of the MPs 
and therefore a lower number of different MPs leads to 
a lower number of zeros in g(cr)^. This can lead to the 



integration volume being more evenly distributed over the 
sphere, thus yielding a higher amount of entanglement. 

Excluding the Dicke states with their infinite amount 
of CPPs, one observes that highly entangled states tend 
to have a large number of CPPs. The prime example for 
this is the case of five qubits, where the classical solution 
with only three CPPs is less entangled than the 'square 
pyramid' state that has five CPPs. In Theorem [T] it was 
shown that 2n — 4 is an upper bound on the number of 
CPPs of positive symmetric n qubit states. For all of 
our numerically determined maximally entangled states - 
including the non-positive ones - this bound is obeyed, 
and for most states the number of CPPs is close to the 
bound {n = 5,8) or coincides with it (n — 4,6,7,12). 
This raises the question whether this bound also holds 
for general symmetric states. To date, neither proof nor 
counterexample is known. 

When viewing the n MPs of a symmetric state as the 
edges of a polyhedron, Euler's formula for convex poly- 
hedra yields the upper bound 2n — 4 on the number of 
faces. This bound is strict if no pair of MPs coincides 
and all polyhedral faces are triangles. Intriguingly, this 
bound is the same as the one for CPPs mentioned in 
the previous paragraph, and the polyhedral faces of our 
numerical solutions come close to the bound (n — 5, 8, 11) 
or coincide with it (n — 4, 6, 7, 10, 12). The faces of the 
polyhedron associated with the Majorana representation 
might therefore hold the key to a proof for 2n — 4 being 
the upper bound on the number of CPPs for all symmetric 
states (with only the Dicke states excluded). 

The case of ten qubits seems to be the first one where 
the maximally entangled symmetric state cannot be cast 
as a positive state. For n — 10, 11, 12, our candidates for 
maximal entanglement are real states, so the question re- 
mains whether the maximally entangled state can always 
be cast with real coefficients. We consider this unlikely, 
firstly because of the higher amount of MP freedom in 
the general case and secondly because many of the so- 
lutions to the classical problems cannot be cast as real 
MP distributions for higher n. For Thomson's problem, 
the first distribution without any symmetry (and thus no 
representation as a real state) arises at n = 13, and for 
Toth's problem at n = 15. 

Upper and lower bounds on the maximal entanglement 
of symmetric states have already been discussed in Section 
with a new proof for the upper bound given in Ap- 
pendix [X] Stronger lower bounds can be computed from 
the known solutions to Toth's and Thomson's problems by 
translating their point distribution into the corresponding 
symmetric state and determining its entanglement. The 
diagram in Figure [14] displays the entanglement of our 
numerical solutions, together with all bounds. 

For n > 5 qubits, the maximally entangled state cannot 
be symmetric or turned into one by LOCC, because the 
lower bound on general states is higher than the upper 
bound of symmetric states. For n — 3, the maximally 
entangled state (W state) is demonstrably symmetric, 
but for n = 4, 5 the numerical solutions for symmetric 
states have less entanglement than the lower bound of 
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FIG. 14. (color online) Scaling of maximal symmetric entan- 
glement with the number of qubits n. The upper bound is 
represented by a black line, while the most entangled Dicke 
states form the lower bound. Their Stirling approximation is 
displayed as a grey line. The most entangled positive symmet- 
ric states found are denoted by blue crosses. For n = 10 — 12, 
the best candidates cannot be cast with positive coefficients 
and are depicted as red stars. The solutions of Toth's and 
Thomson's problems readily provide lower bounds, displayed 
as dashed green lines. The solutions of Thomson's problem are 
generally more highly entangled than those of Toth's problem. 

general states. This would imply that n qubit maximally 
entangled states can be symmetric if, and only if, n < 3. 

VIII. CONCLUSION 

In this paper, we have investigated the maximally entan- 
gled state among symmetric quantum states of n qubits. 
By visualizing symmetric states through the Majorana 
representation and with the help of analytical and nu- 
merical methods, strong candidates for the maximally 
entangled symmetric states of up to 12 qubits were found. 
A comparison with the extremal distributions of Toth's 
and Thomson's problems shows that, in some cases, the 
optimal solution to Majorana's problem coincides with 
that of the two classical problems, but in other cases it 
significantly differs. 

Lower and upper bounds show that the maximal entan- 
glement of permutation-symmetric qubit states scales be- 
tween ©(log -/n) and ©(logn) with the number of qubits 



n. With respect to MBQC, these results indicate that, 
although permutation-symmetric states may not be good 
resources for deterministic MBQC, they may be good for 
stochastic approximate MBQC [HI IS] ■ It also gives bounds 
on how much information can be discriminated locally, 
for which explicit protocols are known in some cases (in 
particular for Dicke states) [TU] . 

We remark that, due to the close relationship of 
distance- like entanglement measures [S], the results for 
the geometric measure give bounds to the robustness of 
entanglement and the relative entropy of entanglement, 
which can be shown to be tight in certain cases of high 
symmetry (TUIUS]. 

We finally note that a similar study has been carried 
out, although from a different perspective, in search of 
the 'least classical' state of a single spin-j system 
(which they call the 'queens of quantumness'). There the 
Majorana representation is used to display spin-j states, 
through the well-known isomorphism between a single 
spin-j system and the symmetric state of n = 2j spin- 
1/2 systems. In this context, the most 'classical' state 
is the spin coherent state, which corresponds exactly to 
a symmetric product state in our case (i.e. n coinciding 
MPs). The problem of [5^ is similar to ours in that they 
look for the state 'furthest away' from the set of spin 
coherent states. However, different distance functions are 
used, so the optimization problem is subtly different and 
again yields different solutions. It is in any case interesting 
to note that our results also have interpretations in this 
context and vice versa. 
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Appendix A: Upper bound on symmetric 
entanglement 

A symmetric n qubit state can be written as 

n 

|V>)=^afee-''|5fc) , 

fc=0 

with Cfc G M, afc G [0, 2n) and the normalization condition 
Tlik'^'k ~ 1- Writing the closest product state as |A) = 
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|cr)®" with |cr) = ce|0) +e''^se|l), we obtain 

n 
fe=0 



(Al) 



Using the set of qubit unit vectors Hi and the uni- 
form measure over the unit sphere dB, the squared norm 



of Equation (All can be integrated over the Majorana 
sphere: 



27r 7r 



\{XH)\'dB 



\{X\ij)\'^smed0dip . (A2) 



k)e-Hi 



Taking into account that J^^ e'"^'^dip = for any integer 
TO 7^ 0, one obtains 



27r -TT 

// 



k=0 



sin 6 d9d(p 



c'^"-'=^sf sin^de 



2^E° 

A:=0 

. 2fn\r{k + l)T{n^k + l 

4^ 1. «^ ^. j rrwT^^ — 

" 1 
^ + 1 



r(n + 2) 
4tt 



k=0 



n + 1 



(A3a) 

(A3b) 

(A3c) 
(A3d) 



The equivalence of Equations ( |A3b ) and ( A3c I follows 
from the different definitions of the Beta function [89j. 
Since the mean value of |(A|'0)P over the Majorana sphere 
is 47r/(n + 1), it follows that G(|V'))^ > l/(n + 1), or 
EciW) — ^og2{n, + 1), for any symmetric n qubit state. 

This result was first shown by R. Renner in his PhD 
thesis j46j , using a similar proof that employs an explicit 
separable decomposition of the identity over symmetric 
subspace. The same proof as ours was independently 
found by J. Martin et al. [55] , 

Appendix B: Proof of theorem [l] 

Class (b): We consider states that have a Z-axis rota- 
tional symmetry with minimal rotational angle Lp = 2it /m^ 
TO G N and 1 < m < n. Figure |15| shows an example 
for m = 5. Due to the rotational Z-axis symmetry and 
the reflective symmetry imposed by Theorem [5| the MPs 
are restricted to specific distribution patters. An arbi- 
trary number of MPs can lie on each of the poles, with 
the remaining MPs equidistantly aligned along horizontal 
circles. The figure shows the two principal types of hori- 
zontal circles that can exist. The upper one is the basic 
type for positive states of five qubits, and the lower one 
a special case where two basic circles are intertwined at 
azimuthal angle ±?9 from the position of the single basic 
circle, respectively. All conceivable horizontal circles of 
MPs can be decomposed into these two principal types. 




FIG. 15. An example of a positive state of 18 qubits with a 
rotational symmetry around the Z-ax;is with angle = 27r/5. 
Two MPs lie on the north pole, one on the south pole, five on 
a single basic circle and ten on two intertwined basic circles. 



According to Equation (14), any CPP |cr) maximizes 
the function nr=i K'^l'^*)!' where the \4>i) are the MPs. 
From Lemma [T] it follows that there must be at least one 
positive I a) = ce|0) -I- Sejl). We first derive the maximum 
number of positive CPPs, and from these the upper bound 
for the total number of CPPs can be immediately obtained 
with the help of Lemma [6] 

For a MP distribution with k MPs on the north pole, I 
MPs on the south pole and the remaining n — k — I MPs 
on horizontal circles, the function to maximize is 

f{e) = (a|o>^-(a|i>'[]/ii(0,)n^2(i9.,e.) , 



where hi{6r) represents the factors contributed by a single 
basic circle with m MPs at inclination Or, and /i2(i?s, Os) 
represents the factors contributed by two basic circles 
with 2m MPs intertwined at azimuthal angles zL-dg, and 
inclination 9^. It is easy to verify that 

h2{^., 0s) = c^"c2™ + 2 cos(to,?,)c^s^c^s^ + s2™s^:" . 
From this it is clear that / can be written in the form 



p p 



f{e) = CgSg 2^ a,Cg 



i=0 



i=0 



where the are positive- valued coefficients, and p is the 
number of basic circles (fc -I- / + pm = n). The number of 
zeros of f'{9) in 9 £ (0, tt) gives a bound on the number of 
positive CPPs. The form of f'{9) is qualitatively different 
for m = 2 and to > 2. With the substitution x = tan(^^/2) 
the equation f'{9) = for to = 2 becomes 



aol + hx^'j - apkx'^P+^ = , with 

b, = a,{l + 2i) - a,_i(fc + 2{p -i) + 2) eR . 

This is a real polynomial in x, with the first and last 
coefficient vanishing if no MPs exist on the south pole 
{I = 0) and north pole {k — 0), respectively. Descartes' 
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rule of signs states that the number of positive roots of a 
real polynomial is at most the number of sign differences 
between consecutive nonzero coefficients, ordered by de- 
scending variable exponent. From this and the fact that 
the codomain of x is M"*", we obtain the result that for 
771 = 2 there are at most p — I, p or p + I extrema of f[9) 
lying in 6 £ (0, tt), depending on whether k and I are zero 
or not. 

For TO > 2, we obtain the analogous result 

with Ci = ai^i[k + (p — i)m + to) G , 
and di — ai{l + im) G . 

From Descartes' rule of signs, we find that there exist 
2p— 1, 2p or 2p+l extrema of f{0) in 6* G (0, tt), depending 
on whether k and / are zero or not. 

With these results it is easy to determine the maximum 
number of global maxima of which are the positive 
CPPs. Case differentiations have to be done with regard 
to TO = 2 or to > 2, whether k and / are zero or not and 
whether p is even or odd. Due to the rotational Z-axis 
symmetry, the non-positive CPPs can be immediately 
obtained. For any positive CPP not lying on a pole, there 
are m — 1 other, non-positive CPPs lying at the same 
inclination (cf. Lemma [6]). For to = 2, the maximum 
possible number of CPPs is (n/2)-|-l {n even) or (n-|-l)/2 
(n odd). This is significantly less than that in the general 
case TO > 2, where the maximum number of CPPs is 
2n — 4. Interestingly, this bound decreases to n if at least 
one of the two poles is free of MPs. 

Class (c): All MPs of a positive state must either lie 
on the positive meridian or form complex conjugate pairs 
(cf . Lemma [H]) . From this the optimization function is 
easily derived as 

n 

/(0)=5]a,cr^s^ , 

i=0 



with real a^. Calculating f'{9) yields the condition for 
the extrema: 

with bi — ai+i{i + 1) — ai_i(n — i + 1) . 

From this, the maximum number of CPPs can be derived 
with Descartes' rule. All CPPs are now restricted to the 
positive meridian and the poles, yielding at most (n-|-3)/2 
CPPs for odd n and (n -|- 2)/2 for even n. □ 



Appendix C: Table of entanglement values 



TABLE I. The table lists the known (n = 2, 3) and numerically 
determined (n > 3) values of the maximal entanglement of 
symmetric n qubit states. The left column lists the extremal 
entanglement among positive symmetric states, and, where 
more entangled non-positive symmetric states are known, they 
are displayed in the right column. 



n 




Eg 


2 


1 




3 


21og23-2« 1.170 




4 


log2 3 ^ 1.585 




5 


1.742 268 948(2] 




6 


21og23- 1 « 2.170 




7 


^ 2.298 691396(2] 




8 


^ 2.445 210 159 




9 


^ 2.553 960 277(2] 




10 


^ 2.679 763 092 


« 2.737 432 003 


11 


^ 2.773 622 669 


^ 2.817 698 505 


12 


^ 2.993 524 700 


log2(243/28) ^ 3.117 



The analytic form is known, but is of a complicated form. 



[1] R. Horodecki, P. Horodecki, M. Horodecki, and K. 

Horodecki, Rev. Mod. Phys. 81, 865 (2009). 
[2] W. Diir, G. Vidal, and J. I. Cirac, Phys. Rev. A 62, 

062314 (2000). 
[3] A. Shimony, Ann. NY. Acad. Sci. 755, 675 (1995). 
[4] T.-C. Wei, M. Ericsson, R M. Goldbart, and W. J. Munro, 

Quant. Inf. Comp. 4, 252 (2004). 
[5] M. Hayashi, D. Markham, M. Murao, M. Owari, and 

S. Virmani, Phys. Rev. Lett. 96, 040501 (2006). 
[6] R. F. Werner and A. S. Holevo, J. Math. Phys. 43, 4353 

(2002). 

[7] D. Gross, S. T. Flammia, and J. Eisert, Phys. Rev. Lett. 

102, 190501 (2009). 
[8] M. Van den Nest, W. Diir, A. Miyake, and H. J. Briegel, 

New J. Phys. 9, 204 (2007). 
[9] C. E. Mora et ai, Phys. Rev. A 81, 042315 (2010). 
[10] M. Hayashi, D. Markham, M. Murao, M. Owari, and 



S. Virmani, Phys. Rev. A 77, 012104 (2008). 
[11] D. Cavalcanti, Phys. Rev. A 73, 044302 (2006). 
[12] L. De Lathauwer, B. De Moor, and J. Vandewalle, 

SIAM J. Matrix Anal. Appl. 21, 1324 (2000). 
[13] T. Zhang and G. H. Golub, SIAM J. Matrix Anal. Appl. 

23, 534 (2001). 

[14] E. Kofidis and P. A. Regaha, SIAM J. Matrix Anal. Appl. 
23, 863 (2002). 

[15] H. Wang and N. Ahuja, IEEE Int. Conf. on Pattern 
Recognition, ICPR (Cambridge, UK) pp 44-47 (2004). 

[16] G. Ni and Y. Wang, Math. Comput. ModeUing 46, 1345 
(2007). 

[17] V. De Silva and L. H. Lim, SIAM J. Matrix Anal. Appl. 

30, 1084 (2008). 
[18] E. D'Hondt and P. Panangaden, Quant. Inf. Comp. 6, 

173 (2006). 

[19] D. M. Greenberger, M. A. Horne, A. Shimony, and 



18 



D. 

9, 
E. 
H. 



A. Zeilinger, Am. J. Phys. 58, 1131 (1990). 

R. H. Dicke, Phys. Rev. 93, 99 (1954). 

R. Prevedel et al, Phys. Rev. Lett. 103, 020503 (2009). 

W. Wieczorek ef al, Phys. Rev Lett. 103, 020504 (2009). 

D. Markham, |arXiv: 1001.03431 (2010). 

G. Toth and O. Giihne, Phys. Rev. Lett. 102, 170503 

(2009) . 

R. Hiibener, M. Kleinmann, T.-C. Wei, C. Gonzalez- 
Guillen, and O. Giihne, Phys. Rev. A 80, 032324 (2009). 
T. Bastin et al, Phys. Rev. Lett. 103, 070503 (2009). 
P. Mathonet et al, Phys. Rev. A 81, 052315 (2010). 

E. Majorana, Nuovo Cimento 9, 43 (1932). 

V. Vedral and M. B. Plenio, Phys. Rev. A 57, 1619 
(1998). 

T.-C. Wei and P. M. Goldbart, Phys. Rev. A 68, 042307 
(2003). 

H. Barnum and N. Linden, J. Phys. A; Math. Gen. 34, 
6787 (2001). 

D. C. Brody and L. P. Hughston, J. Geom. Phys. 38, 19 
(2001). 

R. F. Werner and M. M. Wolf, Phys. Rev. A 64, 032112 
(2001). 

G. Vidal and R. Tarrach, Phys. Rev. A 59, 141 (1999). 
Markham, A. Miyake, and S. Virmani, New J. Phys. 
194 (2007). 

Jung et al, Phys. Rev. A 77, 062317 (2008). 
Zhu, L. Chen, and M. Hayashi, New J. Phys. 12, 
083002 (2010). 

M. B. Plenio and S. Virmani, Quant. Inf. Comp. 7, 1 
(2007). 

G. Toth, J. Opt. Soc. Am. B 24, 275 (2007). 

R. Oriis, S. Dusuel, and J. Vidal, Phys. Rev. Lett. 101, 
025701 (2008). 

J. K. Korbicz, J. I. Cirac, and M. Lewenstein, Phys. Rev. 
Lett. 95, 120502 (2005). 

J. K. Korbicz et al, Phys. Rev. A 74, 052319 (2006). 
J. K. Stockton, J. M. Geremia, A. C. Doherty, and 

H. Mabuchi, Phys. Rev. A 67, 022112 (2003). 

M. Hayashi, D. Markham, M. Murao, M. Owari, and 

S. Virmani, J. Math. Phys. 50, 122104 (2009). 

T.-C. Wei and S. Severini, J. Math. Phys. 51, 092203 

(2010) . 

R. Renner, ETH Zurich, PhD Thesis, l arXiviquanF] 
_ph/0512258 (2005). 
A trivial example of an n qubit state with Eg = (n/2) 
are (n/2) bipartite Bell states, each of which contributes 1 
ebit. Another example is the 2D cluster state of n qubits, 
which has Eq = (n/2) [35]. 
H. Bacry, J. Math. Phys. 15, 1686 (1974). 
J. Zimba and R. Penrose, Stud. Hist. Phil. Sci. 24, 697 
(1993). 

R. Penrose and W. Rindler, Spinors and Space-Ttme. 
Vol.1: Two-Spinor Calculus and Relativistic Fields (Cam- 
bridge University Press, Cambridge, 1984). 
J. Zimba, Electron. J. Theor. Phys. 3, 143 (2006). 
O. Giraud, P. A. Braun, and D. Braun, New J. Phys. 12, 
063005 (2010). 

D. Markham and V. Vedral, Phys. Rev. A 67, 042113 
(2003). 

J. H. Hannay, J. Phys. A: Math. Gen. 29, LlOl (1996). 
J. H. Hannay, J. Phys. A: Math. Gen. 31, L53 (1998). 
P. Leboeuf, J. Phys. A: Math. Gen. 24, 4575 (1991). 



[57; 

[58 

[59 

[60 

[61 

[62 

[63 

[64 

[65 

[66 
[67 

[68 

[69 
[70' 
[71 
[72 
[73' 
[74 

[75 

[76 

[77 
[78' 
[79; 

[80 

[81 



[82' 
[83 



[85 

[86 
[87 



P. Ribeiro, J. Vidal, and R. Mosseri, Phys. Rev. E 78, 
021106 (2008). 

R. Barnett, A. Turner, and E. Demler, Phys. Rev. Lett. 
97, 180412 (2006). 

R. Barnett, A. Turner, and E. Demler, Phys. Rev. A 76, 
013605 (2007). 

R. Barnett, S. Mukerjee, and J. E. Moore, Phys. Rev. 
Lett. 100, 240405 (2008). 

H. Makela and K.-A. Suominen, Phys. Rev. Lett. 99, 
190408 (2007). 

P. Kolenderski and R. Demkowicz-Dobrzanski, Phys. 
Rev. A 78, 052333 (2008). 

P. Kolenderski, Open Systems & Information Dynamics 
17(2), 107 (2009). 

S. Tamaryan, T.-C. Wei, and D. Park, Phys. Rev. A 80, 
052315 (2009). 

L. Chen, A. Xu, and H. Zhu, Phys. Rev. A 82, 032301 
(2010). 

L. L. Whyte, Am. Math. Mon. 59, 606 (1952). 

P. M. L. Tammes, Recueil des trav. bot. neerlandais 27, 

1 (1930). 

E. Weisstein, MathWorld 'Spherical Code', http:// 

mathworld. wolf ram. com/SphericalCode .html 

J. J. Thomson, Phil. Mag. 7, 237 (1904). 

P. Leiderer, Z. Phys. B 98, 303 (1995). 

E. J. Davis, Aerosol Sci. Technol. 26, 212 (1997). 

C. J. Marzec and L. A. Day, Biophys. J. 65, 2559 (1993). 
A. D. Dinsmore et al. Science 298, 1006 (2002). 

H. W. Kroto, J. R. Heath, S. C. O'Brien, R. F. Curl, and 
R. E. Smalley, Nature 318, 162 (1985). 
M. J. W. Dodgson and M. A. Moore, Phys. Rev. B 55, 
3816 (1997). 

T. Erber and G. M. Hockney, J. Phys. A: Math. Gen. 24, 
L1369 (1991). 

J. Leech, Math. Gazette 41, 81 (1957). 

N. Ashby and W. E. Brittin, Am. J. Phys. 54, 776 (1986). 

E. L. Altschuler, T. J. WiUiams, E. R. Ratner, F. Dowla, 

and F. Wooten, Phys. Rev. Lett. 72, 2671 (1994). 

N. J. A. Sloane et al, 'Spherical Codes', http:// 

www2 .research. att . com/ "nj as/packings/ 

D. J. Wales and S. Ulker, 'Global Minima for the Thom- 
son Problem', online at http://www-wales.ch.cam.ac. 
uk/^wales/CCD/Thomson/table .html 

A. A. Berezin, Am. J. Phys. 53, 1036 (1985). 

C. S. Ogilvy and L. Moser, Am. Math. Mon. 58, 492 

(1951). 

ao can be set positive by means of the global phase, and 
for symmetric states with |0) as a GPP it is easy to 
verify that ai = is a necessary condition for the partial 
derivatives of being zero at \uj) — jO). 

The ceiling function \x] is the smallest integer not less 
than X. 

E. Marx, J. Franklin Inst. 290, 71 (1970). 

K. Schiitte and B. L. Van der Waerden, Math. Ann. 123, 
96 (1951). 

J. Martin, O. Giraud, P. A. Braun, D. Braun, and 
T. Bastin, Phys. Rev. A 81, 062347 (2010). 
M. Abramowitz and I. A. Stegun, Pocketbook of Math- 
ematical Functions (Verlag Harri Deutsch, Frankfurt, 
1984). 



